Vector analysis on fractals and applications 



Michael Hinz 1 and Alexander Teplyaev 2 

Abstract. The paper surveys some recent results concerning vector analysis on fractals. 
We start with a local regular Dirichlet form and use the framework of 1-forms and derivations 
introduced by Cipriani and Sauvageot to set up some elements of a related vector analysis 
in weak and non-local formulation. This allows to study various scalar and vector valued 
linear and non-linear partial differential equations on fractals that had not been accessible 
before. Subsequently a stronger (localized, pointwise or fiberwise) version of this vector 
analysis can be developed, which is related to previous work of Kusuoka, Kigami, Eberle, 
Strichartz, Hino, Ionescu, Rogers, Rockner, and the authors. 
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1. Introduction 

In the present article we survey some recent results concerning vector analysis based 
on symmetric local regular Dirichlet forms on locally compact separable metric spaces, cf. 
|16) 1291 I70L 176] . The notions and results we discuss have been introduced in the papers 
|41) 142), 143] . They are based on the approach to differential 1-forms as proposed by Cipriani 
and Sauvageot in |21] in much greater generality, and later investigated by several authors, 
|20l 122) I40L 147] . The constructions are sufficiently robust to apply to symmetric diffusions 
on fractals such as p.c.f. self-similar sets |52L I53L 182] . nested fractals [69], finitely ramified 
fractals |78l 187] . generalized Sierpinski carpets [SI I65L IIP] , spaces of Berlow-Evans-Laakso 
type |79L 180] , and some random fractals |33L 134] . As they are based on Dirichlet forms they 
also apply to classical situations such as Euclidean spaces, domains with sufficiently regular 
boundary and smooth compact Riemannian manifolds. In these cases we recover well-known 
results. 
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2 HINZ AND TEPLYAEV 

A general theme motivating our studies consists of the questions for which elements of 
differential geometry and vector analysis one can find analogs built solely upon the notion 
of energy and how these analogs can be used to formulate and study physical models on 
non-smooth spaces. The space % of 1-forms as constructed in [21L 122] is a Hilbert space. 
Therefore one can identify 1-forms and vector fields, and furthermore introduce other notions 
of vector analysis, as recently done in |41] (which generalizes earlier approaches to vector 
analysis on fractals, see j52L I55L I64L 1721 1811 187]). This is a part of a comprehensive 
program to introduce and study vector equations on general non-smooth spaces which carry 
a diffusion process (or, equivalently, a local regular Dirichlet form). 

Much of the existing literature on analysis on fractals has been concerned with the 
primary problems of construction diffusions on fractals, [6], I10L I13L 131 1. 1371 1521 161L I62L 1631 
I65L I69L I78[ an d references therein], studying their heat kernel decay, [7), I12L I32L 1491 I5CH 
1561 I60j an d references therein], their potential theory, [H], 1111 H4L 1461 171j and references 
therein] , their spectral properties, P^[2l[31l4l[2^ [3l?l[33l[3^[5^ [5^|4^[6^[6^ [8^^ 
and references therein] and some related elliptic and parabolic partial differential equations, 
[25|. I26L 1441 1451 and references therein]. For some recent physics applications of analysis on 
fractals see pQ, [2], I17L I27L I51L I75L 1831 and references therein] , and for analysis on fractals 
in general see [5], 1531 182j . 

Once a diffusion is known to exist, we may regard its infinitesimal generator as the 
Laplacian A, and employ general functional analytic tools (such as semigroup theory or 
variational methods, |24j ) to solve equations of type Au = f and % = An + F(u), and 
even such of form A$(u) = / or ^ = A$(u), with possibly nonlinear transformations F 
and <3>. Note that these equations do not include analogs of first order operators (gradients). 
However, we would like to investigate scalar equations of type 

(1) div(a(Vu)) = / 
or 

(2) Au + b(Vu) = f 

with possibly nonlinear a and b, or vector equations like for instance the Navier-Stokes system 



(3) 



% + (u ■ V)w - Au + Vp = 0, 



dt 

div u = 



or the magnetic Schrodinger equation 

Pint 

(4) i^- = (-iV - Afu + Vu. 

Previous constructions [52, 63], 1811 185] . of first order operators related to diffusions on 
fractals were rather based on purely probabilistic and point-wise approaches, and perhaps 
for this reason not quite flexible enough to fit into a setup that allows to investigate partial 
differential equations containing first order terms. The machinery of [21], 122] . together with 
further developments in [41L I42L I43L 147] , provides a functional analytic definition of a first 
order derivation (respectively gradient) and a framework suitable for a comfortable analysis 
of problems like (|T[) - (SJ) on fractals. 

It is our aim in this paper to highlight elements of this toolkit and to announce some 
related results. We proceed as follows. In the next section we state our main hypotheses 
and collect some useful facts on Dirichlet forms and energy measures. In Section [3] we review 



VECTOR ANALYSIS ON FRACTALS AND APPLICATIONS 



3 



the basic setup of \21\ [22] and discuss related notions of vector analysis proposed in [41J. 
First applications to scalar valued partial differential equations of types ([1]) and ([2]) are then 
presented in SectionHJ and some results on analogs of (j3J) in Section[5j In Section[6]we discuss 
an apporach to (j3J. We also present the definition of related Dirac operators proposed in 
abstract form in |21] and in pointwise form in[43j. 

2. Dirichlet forms and energy measures 

Let X be a locally compact separable metric space and m a Radon measure on X such 
that each nonempty open set is charged positively. We assume that (SjJ 7 ) is a symmetric 
local regular Dirichlet form on L%(X, m) with core C := J 7 C\Co{X). Endowed with the norm 
||/|| c := £(/) 1//2 + sup x \ f\ the space C becomes an algebra and in particular, 

(5) £(fg) 1/2 <\\f\\c\\9\\c, f,gec, 

sec [16|. For any g, h G C we can define a finite signed Radon measure Y(g.h) on X such 
that 

l\ fdT(g,h)=£(fg,h)+£(fh,g)-£(ghJ), feC, 
Jx 

the mutual energy measure of g and h. By approximation we can also define the mutual 
energy measure T(g,h) for general g,h G J 7 . Note that T is symmetric and bilinear, and 
r(<?) > 0, g G T . For details we refer the reader to [29] . We provide some examples. 



Examples 2.1. 

(i) Dirichlet forms on Euclidean domains. Let X = Q be a bounded domain in R™ with 
smooth boundary dQ and 

£(f,g)= [ VfVgdx, f, g ec°°(n). 

Jn 

If Hq(Q) denotes the closure of C°°(Q) with respect to the scalar product £i(f, g) : = 
£, 2 (n)' t nen ^q(^)) ^ s a l° ca l regular Dirichlet form on L 2 (f2). The 
mutual energy measure of /, g G H^lfl) is given by VfVgdx. 

(ii) Dirichlet forms on Riemannian manifolds. Let X = M be a smooth compact 
Riemannian manifold and 

£(f,g)= [ (df,dg) T * M dvol, f,geC°°{M). 

J M 

Here dvol denotes the Riemannian volume measure. Similarly as in (i) the closure of 
£ in L 2 (M,dvol) yields a local regular Dirichlet form. The mutual energy measure 
of two energy finite functions /, g is given by (df, dg) T , M dvol. 

(iii) Dirichlet forms induced by resistance forms on fractals. Let X be a set and (£, J 7 ) a 
local resistance form on it such that X, endowed with the corresponding resistance 
metric R, is complete, separable and locally compact. For any Borel regular measure 
m on (X, R) such that < m(B(x, r)) < 00, the space (J-nL 2 (X, m),£i) is Hilbert, 
and denoting by T the closure of C (X) n T in it, we obtain a local regular Dirichlet 
form (£, J 7 ) on L 2 (X, m) (see for instance [57j Section 9]). Here we have again used 
the standard notation £ x (f,g) = £(f,g) + (f,g) L2 i X m)- 
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Remark 2.1. In Examples 12.11 (i) and (ii) the energy measures have been absolutely 
continuous with respect to the given reference measure. For diffusions on self-similar fractals 
this is typically not true if we choose the corresponding self-similar Hausdorff type measure as 
reference measure, see for instance [15] or [39] . We may, however, use Kusuoka type measures 
as reference measures to produce absolute continuity, see for instance [411 1491 1551 163L 87J. 

3. 1-forms and vector fields 

Following [211 122] we consider C <S> Bb(X), where £>&(X) denotes the space of bounded 
Borel functions on X. We endow this tensor product with the symmetric bilinear form 



(6) (a (g) b, c (g> d) H : = / bddT(a,c), 

a®b,c® d G C ® Bb(X), let ||-||^ denote the associated seminorm on C ® Bb(X) and write 



ker \\-\\ H 



j^ai ® 6i G C®B b (X) : 



22 a i ® & < 



= o 

H 

(with finite linear combinations). To the Hilbert space % obtained as the completion of 
C ® Bb(X)/ker \\-\\ n with respect to ||-|L we refer as the space of differential 1-forms on X, 
cf. [2111221,1501 157]. 

The space % becomes a bimodule if we declare the algebras C and Bb{X) to act on it as 
follows: For a ® b G C ® B b (X), c G C and d G B b (X) set 

(7) c(a ® b) := (ca) g> b - c ® (ab) 
and 

(8) (a®b)d := a® (bd). 

In [21] and [47] it has been shown that © and (|HJ) extend to well defined left and right 
actions of the algebras C and Bb(X) on H. From ([6]) and the Leibniz rule for energy measures, 
see |29l Theorem 3.2.2], it can be seen that left and right multiplication agree for any c G C, 
and as 

max{||(a (g) 6)c||^ , \\c(a <S> b)\\ n } < sup \c\ \\a <S> b\\ n , 

x 



it follows by approximation that they agree for all c G Bb(X), see 
A derivation operator d : C — > H can be defined by setting 

df:=f ®1. 

It obeys the Leibniz rule, 

(9) d(fg) = fdg + gdf, f,geC, 
and is a bounded linear operator satisfying 

(10) \\dff H = S(f), fee. 

On Euclidean domains and on smooth manifolds the operator d coincides with the classical 
exterior derivative (in the sense of L 2 -differential forms). Details can be found in |21L 1221 
HDl fill [4T], 

Being Hilbert, H is self-dual. We therefore regard 1-forms also as vector fields and d as 
the gradient operator. Let C* denote the dual space of C, normed by 

Hl c *=sup{K/)|:/GC,||/|| c <l}. 
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Given /, g G C, consider the functional 

u^d*(gdf)(u):=-(du,gdf) H = - [ g dT(u, f) 

J x 

on C. It defines an element d*(gdf) of C*, to which we refer as the divergence of the vector 
field gdf. 

Lemma 3.1. The divergence operator d* extends continuously to a bounded linear operator 
from TL into C* with \\d*v\\ c * < ||i>|| w , v G Ti. We have 

d*v(u) = — (du, v) u 

for any u G C and any v G "H. 
The Euclidean identity 

div (g grad f) = gAf + V/Va 

has a counterpart in terms of d and d*. Let (A,dom A) denote the infinitesimal L 2 (X,fi)- 
generator of (£, J r ). 

Lemma 3.2. We have 

d*(gdf) = gAf + T(f,g) , 
for any simple vector field gdf , f,g G C, and in particular, Af = d*df for f G C. 

Proofs of these results are given in |41l Section 3]. This distributional perspective can be 
complemented by the following point of view. The operator d, equipped with the domain C, 
may be seen as densely defined unbounded operator 

d: L 2 (X,m) ->H. 

Since (S, J 7 ) is a Dirichlet form, d extends uniquely to a closed linear operator d with domain 
dom d = J 7 . The divergence d* , seen as an operator 

d* : U -> L 2 {X,m), 

will be unbounded, note that in general the inclusions C C L,2(X,m) C C* are proper. As 
usual v G % is said to be a member of dom d* if there exists some v* G L 2 (X, m) such that 
(u, v*) L , x m s = — (du,v) n for all ueC. In this case d*v := v* and 

d * V )L 2 (X,m) = ~ ( 9u i v )h » M G C ' 

i.e. —(9* is the adjoint operator of 9. It is immediate that {Of : / G dom A} C dom d* . As 
9 is densely defined and closed, the domain dom d* of d* is automatically dense in %. 
We say that (E,F) admits a spectral gap if there exists some c > such that 

(11) f (f-fx) 2 dm<cS(f) 

Jx 

for any / G J 7 , where = J x / dm. If (SjJ 7 ) has a spectral gap, then the image 
7m (9 of 9 is a closed subspace of In this case the space % decomposes orthogonally into 
Imd and its complement (Imd) , what implies (Imd) ± = kerd*, and as a consequence we 
observe the following explicit description of domd*. 
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Corollary 3.1. Assume that {S,^) admits a spectral gap, (TO]). Then the domain 
dom d* agrees with 

{v E H : v = df + w : f <E dom A , w G ker d* } . 
For any v = Of + w with f G dom A and w G ker d* we have d*v = Af. 
The proof is short and straightforward, see |43l Corollary 2.2]. 

4. Scalar PDE involving first order terms 

The results of the preceding section may be used to obtain some results on equations of 
type ([T]) and ([2]). We quote from [41[ Section 4]. First consider the quasilinear equation 

(12) d*a{du) = f. 

on L 2 (X,m). In the situation of Example I2.1l (i) it agrees with ([T]). Assume that a : 7-L — > H 
satisfies the following monotonicity, growth and coercivity conditions: 

(13) i a ( v ) ~ a(w),v — w) n > for all v , w G Im d, 

(14) \\a(v)\\ H < cq(1 + \\v\\ H ) for all v e Im d 
with some constant Cq > 0, and 

(15) (a(v),v)n,> Cx\\v\\ n — c 2 for all v G Im d 

with constants c\ > 0, c 2 > 0. For simplicity we assume the validity of a Poincare inequality, 

(16) Wf\\l 2 (x, m) <cpS(f) 

with some constant Cp > for all / G L 2 (X,m) with j x fdm = 0. A function u G J 7 is 
called a weak solution to [W\! if 

(a(du),dv) H = - (f,v) L2{Xm) for all v G T. 

By classical methods, [241 Section 9.1], we obtain the following result. 

Theorem 4.1. Assume a satisfies (T73]) ; fT^] ) and (T73J) and suppose ( li6J) holds. Then 
[W\l has a weak solution. Moreover, if a is strictly monotone, i. e. 



(17) ( a ( w ) ~ a(w), v — w) n > c 3 \\v — /or a// 1> , u> G /m 5 
with some constant C3 > ; then [W\! has a unique weak solution. 

An analog of ([2]) can be treated in a similar manner. Consider 

(18) - Au + 6(0u) + qu = 0, 

where > and 6 is a generally non-linear function- valued mapping on "H. Assume that 
b : "H — > L 2 (X, m) is such that 

(19) H^)llL 2 (x, m) < c 5(l + IMU, velmd, 
with some C5 > 0. A function « € J is called a weak solution to ffTSj) if 

£(w,t>) + (b(du),dv) H + g (w,f)L 2 (x,m) = for a11 v ^ 
From [24[ Section 9.2.2, Example 2], we then obtain the following. 

Theorem 4.2. Assume that the embedding J 7 C L 2 (X,m) is compact and that [W\l holds. 
Then for any sufficiently large g > there exists a weak solution to (T72J]. 
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5. Navier-Stokes equations 

In this section we comment on equations of type ([3]) which provide some more interesting 
applications for the notions discussed in Section [3j 

Together with suitable boundary conditions the Navier-Stokes system ([3]) describes the 
flow of an incompressible and homogeneous fluid in a Euclidean domain with velocity field u 
and subject to the pressure p. In a one-dimensional situation it reduces to an Euler equation 
du/dt+dp/dx = that has only stationary solutions. In |42] we have proposed to investigate 
an analog of ([3]) on compact connected topologically one- dimensional fractals X. We collect 
some items necessary to formulate it. 

Assume that the space X is compact, connected and topologically one- dimensional and 
that (SjJ 7 ) admits a spectral gap, (ITT|) . Combined with several results on Hodge decom- 
positions and topology, cf. [421 Sections 4,5 and 6], the assumption of topological one- 
dimensionality had motivated to define a Laplacian Ai on 1-forms by 

(20) Ai := dd*, 

seen as an unbounded operator on H with domain dom A x = {u G dom d* : d*u G J 7 }. 

Theorem 5.1. The operator (Ai, dom Ai) is a self-adjoint operator onH. 

A proof can be found in [43[ Section 6]. Theorem 15.11 allows to talk about harmonic 
forms: A 1-form u G H, is called harmonic if u G dom Ai and A±u = 0. From compactness 
and topological one-dimensionality can deduce the following, cf. [43[ Theorem 6.2]. 

Theorem 5.2. A 1-form u G "H is harmonic if and only if it is in (Im d)- 1 . 

The proof of of Theorem 15. 21 is rather subtle, it involves a description of (Imd) 1 - in terms 
of locally harmonic forms. We refer the reader to [42J. Note that Theorem 15.21 indicates that 
in this situation the definition (1201) is appropriate. 

Also for the convection term (u ■ V)w in ([3]) we propose a substitute which by one- 
dimensionality seems reasonable. Our choice is motivated by the Euclidean situation: Given 
a vector field u, the quantity 

(21) - J \u\ 2 divvdx, 

seen as a functional on a space of test vector fields v, provides a formulation of V|w| 2 in the 
weak sense. In our situation we set 

dom c d* :={ve dom d* : d*v G C(X)} 

and given u G H, define 

(22) dT H (u)(v) := - ((d*v)u,u) H , v G dom c d*. 

This seems reasonable by a fiberwise (respectively m-a.e. pointwise) representation for 7-L 
proved in |41l Section 2] and |43l Theorem 2.2]: 

Theorem 5.3. Let v be a Radon measure such that all energy measures are absolutely 
continuous with respect to v . There are a family of Hilbert spaces {H x } xeX an< ^ surjective 
linear maps u t— > u x from 7i onto T-i x such that the direct integral 'H. x v{dx) is isometrically 
isomorphic to H and in particular, 

\\ UJ \\u = / W^An x v(dx), a; G T-L. 
Jk 
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Theorem 15.31 itself is more general, it does neither require X to be compact or topologi- 
cally one- dimensional nor {S,^) to admit a spectral gap. 

If we replace the Euclidean norm | • | in (I2T]) by the norms ||-|| Wx . of the fibers H x , we 
arrive at ( I22p . On the other hand, we have the classical identity 

-V|w| 2 = (u ■ V)m + u x curlw. 

In a one- dimensional situation there should be no nonzero 2-forms, hence curlw should be 
trivial, so that ( 122]) is a good substitute for {u ■ V)w. See [42J for a more detailed discussion. 
Altogether this gives a strong heuristic motivation to regard 



(23) 



§ + ±dr n (u)-A lU + d P = o 

d*u = 0. 



as a suitable analog of a (131) on a compact topologically one-dimensional space. Note that 
this is a boundary free formulation. We say that a square integrable dom <9*-valued function 
u on [0, oo) provides a weak solution to with initial condition uq G ker d* if 



(24) 



(u(t),v) H - (uo,v) n + /„* dT n (u(s))(v)ds + jl (d*u(s),d*v} L2(Xm) ds = 
d*u{t) = 

for a.e. t G [0, oo) and all w G fcer 9*. By some immediate simplifications we then observe 
stationarity and uniqueness of solutions. In other words, the behavior of the system on a 
compact topologically one-dimensional space resembles its behavior on a one-dimensional 
Euclidean space. 

Theorem 5.4. Any weak solution u of K2$\) is harmonic and stationary, i.e. u is 
independent of t G [0, oo). Given an initial condition u$ the corresponding weak solution is 
uniquely determined. 

Note that we have not made any restriction on the Hausdorff dimension du of X. Indeed 
there are examples of spaces of any Hausdorff dimension 1 < d H < oo such that the previous 
Theorem holds. It is the topological dimension that governs the behavior of fl2Ul) . 

Remark 5.1. Logically Theorem 15.31 is not needed to set up the model (|23|) . we have 
just included it here to support the intuition behind our choice of substitute terms. We 
would also like to remark that even though the energy measures might not be absolutely 
continuous with respect to the initial reference measure m, one can always construct a finite 
Radon measure v that has this property. 

For the rest of this section we specialize further to the situation of Examples 12. II (iii) . that 
is, we assume X to be a set and (S,^) to be a resistance form on it such that X, together 
with the resistance metric R, is a compact and connected metric space. We further assume 
that m is a Borel regular measure on X as in Examples 12.11 (iii) so that consequently a local 
regular Dirichlet form (S, J 7 ) is obtained by taking the closure. Then all our previous results 
may be applied for (S,^). We finally assume that (X,R) is topologically one- dimensional. 

Remark 5.2. We conjecture that any set that carries a regular resistance form is a 
topologically one-dimensional space when equipped with the associated resistance metric. 

Note that in the resistance form case points have positive capacity. This property allows 
to prove the following equivalence, [42, Section 5]. 
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Theorem 5.5. Assume that points have positive capacity and topological dimension is 
one. Then a nontrivial solution to ( ES| ) exists if and only if the first Cech cohomology i? 1 (X) 
of X is nontrivial. 

In the resistance form context Neumann derivatives are well-defined, and it is not difficult 
to see that if the Navier-Stokes system (I2"3"j) is considered with a nonempty boundary, it may 
have additional nontrivial solutions arising from solutions of a related Neumann problem. 

Let B C X be a finite set, which is interpreted as the boundary of X. By Gb we denote 
the Green operator associated with the boundary B with respect to (£, J 7 ), |54[ Definition 
5.6], and T>^ its image in T. Let 'Kb denote the 5-harmonic functions with respect to 
(£, J 7 ), |54[ Definition 2.16], and note that T = J ' B © 'Kb-, where 

Tb '■= {u G T : u\b = 0} . 

A 5-harmonic function h is harmonic on B c in the Dirichlet form sense, more precisely, it 
satisfies S(h, ip) — for all ip G J- b- The space V L := V L Ba + H B is seen to be independent of 
the choice of B, [54] Theorem 5.10]. For any u G V L and any p G X the Neumann derivative 
(du) p of u at p can be defined, |54[ Theorems 6.6 and 6.8]. If is a function on B, then a 
function h v G J 7 is called a solution to the Neumann problem on B c with boundary values 
if if it is harmonic on B c and satisfies (dh) p = <p(p) for all p G B. Such a Neumann solution 
h v exists and is unique if and only if ip is such that 

5>(p) = o. 

pes 

We use the notation 7-L(B c ) = closspan {v 1b c '■ v G H}. A square integrable dom d*- 
valued function u on [0, oo) provides a weaA; solution to $E3j) on B c if 

(25) J _ (U(0) ' V) « + ^ 5r «( u ( s ))(^) rfs + Jo ( d * u (s), d*v) L2(X!m) ds = 

for a.e. t G [0, oo), all v G dom d* (1 H(B C ) and all ip G 

Theorem 5.6. Assume that points have positive capacity and topological dimension is 
one. If h is the unique, up to an additive constant, harmonic function on B c with normal 
derivatives <p on B, then 

u(t) = dh, t G [0,oo). 
is the unique weak solution to / pggj) on B c with the Neumann boundary values ip on B. 

Remark 5.3. In (|25|) we have considered weak solutions to ( 123]) . For weak solutions the 
pressure p does not occur explicitely. However, any definition of strong solution to |3] should 
lead to the relation 

p(jt) = ~T(h), tG[0,oo), 

seen as an equality of measures. A similar statement could be written for the boundary free 
case (1241). 

For more details and the Hodge theory leading to the statements of this section we refer 
the reader to |43j . 
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6. Magnetic Schrodinger equations 

We turn to results concerning the magnetic Schrodinger equation fll]). To discuss this 
equation we do not need to assume that X is compact or topologically one- dimensional. As in 
Section [2] it may just be an arbitrary locally compact separable metric space equipped with a 
Radon measure m that charges any nonempty open set positively and carrying a symmetric 
local regular Dirichlet form (£, T) on L 2 (X,m) with core C := J 7 fl C (X). However, to 
investigate (J4]) we will now assume that (SyJ 7 ) possesses energy densities with respect to the 
reference measure m, i.e. for any g,h G J 7 the measure T(g, h) is absolutely continuous with 
respect to m. 

Remark 6.1. As previously mentioned in Remark 15.1 1 we can always construct a measure 
v with respect to which all energy measures are absolutely continuous. For the cases that 
the given Dirichlet form (£, J 7 ) on L2(X,m) is transient or induced by a resistance form we 
have shown in |41] that (£, C) is a closable form on the space L 2 (X, v) of functions that are 
square integrable with respect to this new measure v. Then this change of measure merely 
amounts to a change of domains. It is not difficult to show that under mild assumptions 
(£, C) is always closable with respect to this measure v. We will discuss this matter in a 
later paper. 

In |43j we have studied analogs of the magnetic Hamiltonian (— zV — A) 2 + V and in 
particular, have verified their essential self-adjointness. To sketch this result, let L 2 ^(X 1 m), 
jFc, Cc and %c denote the natural complexifications of L 2 (X,m), J 7 , C and "H, respectively. 
The natural extensions to the complex case of 8 and the corresponding energy measures 
T(f,g) are again denoted by the same symbol. Note that they are conjugate symmetric and 
linear in the first argument. If both arguments agree, they yield a real nonnegative number 
and a real nonnegative measure, respectively. The first result concerns related quadratic 
forms. Here we use the notion of quadratic form in the sense of [73[ Section VIII. 6]. 

Proposition 6.1. Let a eW and V G L^X^m). The form £ a > v , given by 

£ a ' V (f, 9) = ({-id - a)f, (-id - a)g) n + (fV, g) L2{x , m) , f,ge C c , 

defines a quadratic form on L 2j c(A, m). 

Proposition [6J] is a slight variation of Proposition 4.1 in |43] and up to inessential details 
it has the same proof. Here a is seen as the magnetic vector potential replacing A in (J3|) and 
V is the electric scalar potential. 

Now recall the fiberwise representation of H. from Theorem 15.31 We define the space of 
real vector fields of bounded length by 

Hog ■= jt> = (v x ) xeX e H : \\v.\\ Ht . E Loo(A, m)| . 

If the potential a is recruited from Hoo then we can obtain the closedness of S a,v and 
the essential self-adjointness of the associated operator from straightforward perturbation 
arguments, cf. [43[ Theorem 4.1]. Recall that A denotes the generator of (S,^). We denote 
its complexification by the same symbol. 

Theorem 6.1. Let a e and V G L OQ (X,m). 
(i) The quadratic form (£ a,v ,J : ' c ) is closed. 
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(ii) The self-adjoint non-negative definite operator on L 2j c(A, m) uniquely associated 
with (£ a,v , J-'c) is given by 



The operator H a ' is a natural generalization of the quantum mechanical Schrodinger 
Hamiltonian (— zV — A) + V from (j3J). By Theorem 16.11 we have established a suitable 
framework to study a fractal counterpart 




of the evolution equation ((4]). 

As they are closely related to magnetic Hamiltonians, we conclude this section by a brief 
look at Dirac operators. We have introduced a local Dirac operator in [43] . Up to sign and 
complexity conventions it is defined as a matrix operator 



acting on H © Hi = L 2 (X, m) © "H. We consider D as an unbounded linear operator with 
domain dom D := J 7 © dom d* and have the following result, obtained in abstract form in 
|21j . and in pointwise form in [43]. 

Theorem 6.2. The operator (D, dom D) is self-adjoint operator on H © Hi. 

Note that as a consequence we also obtain a local matrix Laplacian D 2 acting on Hq(BHi. 

According to |42L 143] . this Dirac operator is naturally related to the topological structure 
of the fractals space and, in a certain natural a sense, to the differential geometry of the 
fractal (see |18L 1191 and references therein] for a discussion of the notion of a Dirac operator 
in the context of non- commutative analysis). In particular, for topologically one dimensional 
fractals (of arbitrary Hausdorff and spectral dimensions) our Dirac operator gives rise to a 
natural Hodge Laplacian dd* + d*d on the appropriate differential complex. It will be 
the subject of future work to study the Hodge Laplacian for higher order differential forms 
defined in the probabilistic or Dirichlet form sense. 

As a side remark we note that sometimes there may be a different convention for the 
Dirac operator in a complex setup. For instance 



has signs and imaginary factors which are somewhat more suitable in relation to the magnetic 
magnetic Schrodinger operator H a,v . 

Finally, we would like to point out related perturbation results. Assume that b E H^ 
and set 



/, g G J 7 . Here b x and d x f denote the images of b and df under the projection from H onto 
H x as in Theorem 15.31 For a > write 



H a,v = (_ id _ a y(_ id - a ) + V, 
and the domain of the operator A is a domain of essential self- adjointness for H a ' 



(26) 






Q(f,g) ■= Q{f,9) + a{f,g) La{x , m) - 



We may then conclude the following. 
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Theorem 6.3. 

(i) For any a > the form (Q a , J 7 ) is closed on L2(X,m). It generates a strongly 
continuous semigroup of bounded operators on /^(X, m). 

(ii) If a > is sufficiently large then the associated semigroup is positivity preserving. 

(iii) The generator L Q of Q is given by 

L Q u(x) = Au(x) + (b x , d x f) n , u 6 dom A. 
See |41[ Section 10] and the references therein, in particular |28j. 
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